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ABSTRACT 

Using self-consistent three-dimensional (3D) IV-body simulations, we investigate the physical prop¬ 
erties of non-axisymmetric features in a disk galaxy created by a tidal interaction with its companion. 
The primary galaxy consists of a stellar disk, a bugle, and a live halo, corresponding to Milky-Way 
type galaxies, while the companion is represented by a halo alone. We vary the companion mass and 
the pericenter distance to explore situations with differing tidal strength parameterized by either the 
relative tidal force P or the relative imparted momentum S. We find that the formation of a tidal 
tail in the outer parts requires P > 0.05 or S > 0.07. A stronger interaction results in a stronger, 
less wound tail that forms earlier. Similarly, a stronger tidal forcing produces stronger, more loosely 
wound spiral arms in the inner parts. The arms are approximately logarithmic in shape, with both 
amplitude and pitch angle decaying with time. The derived pattern speed decreases with radius and 
is close to the Q — n/2 curve at late time, with D and n denoting the angular and epicycle frequencies, 
respectively. This suggests that the tidally-induced spiral arms are most likely kinematic density waves 
weakly modified by self-gravity. Compared to the razor-thin counterparts, arms in the 3D models are 
weaker, have a smaller pitch angle, and wind and decay more rapidly. The 3D density structure of 
the arms is well described by the concentrated and sinusoidal models when the arms are in the non¬ 
linear and linear regimes, respectively. We demonstrate that dynamical friction between interacting 
galaxies transfers the orbital angular momentum of one galaxy to the spin angular momentum of the 
companion halo. 

Subject headings: galaxies: spiral — galaxies: structure — galaxies: interactions — galaxies: evolution 
— methods: numerical 


1. INTRODUCTION 

Spiral arms play an important role in galactic evolution 
in diverse ways (e.g., Buta & Combes 1996; Kormendy 
& Kennicutt 2004; Buta 2013; Sellwood 2014 and refer¬ 
ences therein). They exert non-axisymmetric torque to 
stars and gas clouds in galaxy rotation and cause their 
radial migrations, leading to secular density changes of 
the disks (e.g., Foyle et al. 2010; Roskar et al. 2012; Baba 
et al. 2013; Kim & Kim 2014). They also provide sites 
for active star formation, which is either triggered or or¬ 
ganized by the stellar spiral potentials (e.g., Elmegreen 
& Elmegreen 1986; Elmegreen 1995; Sleath & Alexan¬ 
der 1996; Bertin & Lin 1996; Seigar & James 2002; see 
also McKee & Ostriker 2007 and Dobbs & Baba 2014 
for review). 4 Therefore, understanding the nature and 
properties of spiral arms is crucial to understand secu¬ 
lar, chemical, and dynamical evolution of disk galaxies. 

In terms of lifetime of spiral arms, the theory of spi¬ 
ral structure has forked into two branches: long-lived 
quasi-stationary density waves (Lin & Shu 1964, 1966) 
and short-lived transient features (Toomre 1969; Toomre 
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& Toomre 1972). The first picture requires that self¬ 
gravity plays a key role in organizing density waves into 
a self-sustained global pattern that rotates almost rigidly 
about the galaxy center (Bertin et al. 1989a, b; Lowe et al. 
1994; Bertin & Lin 1996). It successfully predicts system¬ 
atic offsets between narrow dust lanes and star-forming 
regions associated with the arms (e.g., Roberts 1969; see 
also Section 6.4.3 of Binney & Tremaine 2008), although 
it does not address the origin of the density waves. In 
the second picture, on the other hand, spiral arms are 
transient waves that wind up over time as they propa¬ 
gate inward in the radial direction. In this case, spiral 
arms last only for a few rotation periods. Numerical 
simulations of isolated disk galaxies show that noises in¬ 
herent in a stellar disk (e.g., Fujii et al. 2011; Grand et 
al. 2012, 2013; Baba et al. 2013) or perturbations pro¬ 
vided by giant molecular clouds (e.g., D’Onghia et al. 
2013) are amplified as they swing from leading to trail¬ 
ing configurations (e.g., Goldreich & Lynden-Bell 1965; 
Julian & Toomre 1966; Toomre 1981), and form transient 
but recurrent spiral arms. Such swing-amplified arms in 
isolated galaxies are usually ragged with multiple arms 
rather than being grand-design spirals with prominent 
two arms. 

Observations indicate that among samples that include 
both non-barred and barred spiral galaxies, the probabil¬ 
ity of having grand-design arms is higher for galaxies with 
companions (e.g., Kormendy & Norman 1979; Elmegreen 
& Elmegreen 1982, 1987). 5 . Recently, Kendall et al. 

5 Having a bar increases the probability to possess grand-design 
arms especially in binary systems (Elmegreen & Elmegreen 1982). 
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(2015) analyzed spiral structures in a sample of galaxies 
from the Spitzer Infrared Nearby Galaxies Survey and 
found that the strength of grand-design arms is rather 
tightly correlated with tidal forcing from nearby com¬ 
panion galaxies, while the arm morphologies depend very 
weakly on the galaxy parameters such as stellar mass, gas 
fraction, disk/bulge ratio, rotational velocity, etc. This 
indicates that a large fraction of grand-design arms, es¬ 
pecially for very strong arms, are driven most likely by 
tidal interactions. 

Toomre & Toomre (1972) pioneered a numerical study 
for tidal interaction of disk galaxies using noninteract¬ 
ing test particles, and found that extended structures 
such as tidal tails and bridges develop in the course 
of galaxy encounters. Subsequent iV-body simulations 
focused on very strong encounters that lead to galaxy 
mergers or significant transformation of disk morpholo¬ 
gies (e.g., Farouki & Shapiro 1982; Hernquist 1990b; 
Hernquist 1992; Barnes 1992; Mihos & Hernquist 1994; 
Barnes 1998; Naab & Burkert 2003; Cox et al. 2006; Chan 
& Junqueira 2014). Some studies concentrated on the or¬ 
bital parameters required to reproduce the observed mor¬ 
phologies and kinematic features of interacting galaxies 
(e.g., Howard & Byrd 1990; Elmegreen et al. 1991; Salo 
& Laurikainen 1993, 2000a), while others investigated 
properties of tidally-induced spiral arms (e.g., Sundelius 
et al. 1987; Byrd & Howard 1992; Donner & Thomasson 
1994; Salo & Laurikainen 2000b). In particular, Sun¬ 
delius et al. (1987) found that the spiral arms generated 
in a cold disk by weak tidal interaction wind up from Sc 
to Sa appearances. Byrd & Howard (1992) showed that 
inner spiral arms are created if the tidal strength param¬ 
eter P (see below for definition) is larger than 0.01, sug¬ 
gesting that even a low-mass companion can tidally ex¬ 
cite grand-design arms if the pericenter distance is small 
enough. To characterize the tidal strength, Elmegreen 
et al. (1991) instead used the S parameter (see also be¬ 
low for definition) that takes into account the interaction 
duration, finding that a tidal encounter with S > 0.019 
deforms the outer disk into an “ocular” shape. 

Numerical simulations often show that tidally-induced 
arms are transient and posses the characteristics of kine¬ 
matic density waves. For example, Donner & Thomas¬ 
son (1994) found that spiral arms driven by tidal forces 
resemble kinematic density waves when self-gravity is 
weak. In modeling tidal interactions of the M51/NGC 
5195 system, Salo & Laurikainen (2000b) found that spi¬ 
ral arms in M51 cannot be described by a single pattern 
speed, indicative of kinematic density waves. Baba et al. 
(2013) showed that spiral arms in an isolated galaxy gen¬ 
erated by swing amplification of random perturbations 
have a pattern speed that decreases with the galactocen- 
tric radius R (see also D’Onghia et al. 2013; Michikoshi & 
Kokubo 2014). Chan & Junqueira (2014) also found that 
the pattern speed of m = 2 arms driven by a satellite on 
a bound eccentric orbit decreases with radius. 

To quantify how the physical properties of the tidally- 
driven arms depend on S , Oh et al. (2008, hereafter Pa¬ 
per I) ran a series of IV-body simulations by considering 
a stellar disk inside a fixed halo interacting with a point 
mass companion on a prescribed parabolic orbit. Paper I 
found that a tidal bridge consists of the disk particles 
pulled out by the tidal perturbations, with their epicy¬ 
cle phases locked to the companion, while a tidal tail 


forms as strongly perturbed near-side particles overtake 
mildly perturbed far-side particles. Paper I also found 
that a stronger encounter produces stronger, more open 
arms that start to develop earlier at smaller R. The arm 
pattern speed turned out to be a decreasing function of 
R even when the arms are strongest, and converges to 
the — k/2 curve as they decay. This suggests that 
tidally-driven arms are unlikely to be quasi-stationary 
density waves. Similar results were obtained by Dobbs 
et al. (2010) who included a gaseous disk as well in sim¬ 
ulating the M51/NGC 5195 system, although gravity of 
the gaseous component tends to increase the arm pat¬ 
tern speed. Struck et al. (2011) suggested that these 
are caustic waves maintained by coherent epicycle oscil¬ 
lations triggered by tidal forcing. 

While the results of the Paper I are informative in 
assessing the quantitative effects of tidal perturbations 
on the arm properties, they were based on highly ideal¬ 
ized galaxy models. First of all, Paper I considered an 
infinitesimally-thin, two-dimensional (2D) stellar disk. 
This not only neglects nonplanar motions but also over¬ 
estimates self-gravity of the disk, making the spiral arms 
stronger than in disks with finite thickness. In addi¬ 
tion, Paper I employed fixed gravitational potentials for 
the halo/bulge and the perturbing companion. This pre¬ 
cludes the possibility of dynamical friction occurring due 
to their gravitational reactions to the tidal perturbations. 
Moreover, the galaxy and the companion in Paper I were 
set to follow the prescribed parabolic orbits, which ig¬ 
nores their orbital decay caused by angular momentum 
loss. By evolving the system in the frame in which the 
galaxy remains stationary, Paper I also ignored indirect 
forces arising from the orbital motion of the galaxy rela¬ 
tive to the center of mass of the system, potentially sup¬ 
pressing the growth of lopsided spiral modes in the stellar 
disk (e.g., Adams et al. 1989; Ostriker et al. 1992). 

To overcome the caveats mentioned above, we in this 
paper extend Paper I by considering a self-consistent 
three-dimensional (3D) galaxy model in which the stel¬ 
lar disk has a finite thickness and a central bugle and a 
dark halo are represented by live particles rather than by 
fixed potentials. The perturbing companion is modeled 
by a live halo alone for simplicity. As in Paper I, we 
control the strength of a tidal encounter by varying two 
parameters: the pericenter distance and the galaxy-to- 
companion mass ratio. Our objective is four-fold. First, 
we wish to measure the physical properties of tidally- 
induced spiral arms in more realistic 3D models and 
compare them with those in the razor-thin counterpart. 
We will also find the quantitative dependence of the arm 
properties on both S and P. Second, Springel & White 
(1999) and Dubinski et al. (1999) showed that the forma¬ 
tion of tidal tails depends on the shape of the galactic po¬ 
tential. We will show that the tail formation depend not 
only on the potential shape but also on the tidal streg- 
nth. Third, we want to explore 3D density structures 
inside spiral arms and compare them with the analytic 
formulae proposed by Cox & Gomez (2002). Finally, we 
will show that dynamical friction of the companion in 
the course of a tidal interaction can be a source of the 
spin angular momentum of the primary halo. 

This paper is organized as follows. In section 2, we 
describe our galaxy model and introduce the parameters 
for tidal strength. In Section 3, we present the proper- 


Tidally-Driven Disk Structures 


3 




Figure 1. (a) Radial profile of the circular velocity of the initial 

disk and the contribution of each component, (b) Radial distribu¬ 
tions of the characteristic angular frequencies fl (solid), C2 — re/2 
(dotted), and f2 + re/2 (dashed). 


ties of tidal tails formed in the outer parts and spiral 
arms in the inner parts. In Section 4, we analyze the 3D 
density structures of spiral arm in the radial and verti¬ 
cal directions and compare them with the prediction of 
Cox & Gomez (2002). In Section 5, we explore the or¬ 
bital decay of the galaxy caused by dynamical friction, 
which results in spin-up of the initially non-rotating halo. 
Finally, we summarize our results and discuss their as¬ 
tronomical implications in Section 6. 


2. MODEL 
2.1. Galaxy Models 

In this paper we use 3D N -body simulations to study 
the tidal interaction of a disk galaxy with its companion. 
The galaxy consists of a stellar disk, a spherical bulge, 
and a dark matter halo; the effect of a gaseous component 
is not considered in the present study. For the halo, we 
adopt a truncated Hernquist (1990a) model 

{ M h r h 

Ph{t)=\ 2tt r(r + rh) 3 ’ “ ^ ( 1 ) 

I 0, for r > rtr, 


with the mass Mh = 5.15 x 10 11 M 0 , the scale radius 
rh = 10 kpc, and the truncation radius ?'t r = 200 kpc, 
within which 90% of Mh is enclosed. For the bulge, we 
take a Plummer sphere 


_ 3Mb rl 

47r (r 2 + r 2 ) 5 / 2 ’ 


( 2 ) 


with Mb = 9.80 x 10 9 M 0 and r b = 0.23 kpc. 


The stellar disk is initially axisymmetric with density 
distribution 


Pd(R, z) 


M d ( R \ 
47 rh 0 R 2 d 6XP \ R d ) 



(3) 


where M d = 5.17 x 10 10 M 0 is the total disk mass, R d = 
3.4 kpc is the radial scale length, and hg = 0.33 kpc is the 
vertical scale height. The corresponding stellar surface 
density is T, d = 37.4 M 0 kpc -2 at R = 10 kpc. The 
mass ratio of each component of the galaxy is 0.9 Mh : 
Mb : M d = 9.1 : 0.2 : 1, with the total mass M ga i = 
5.30 x 10 11 M 0 within the truncation radius. 

Figure 1 plots the radial distribution of the circular 
velocity and the contribution of each component as well 
as the characteristic angular frequencies ft and 0 ± n/2 
in the initial disk. The presence of a strong bulge makes 
the Q — k/2 curve rise steeply toward the center, which 
suppresses the formation of a bar in our tidal encounter 
models (e.g., Toomre 1981; Sellwood 2000). To find the 
equilibrium velocity distribution of the disk particles un¬ 
der the total gravitational potential, we solve the ax¬ 
isymmetric Jeans equations using the method described 
in Hernquist (1993) and Springel et al. (2005). We en¬ 
sure that the dark matter halo and bulge do not rotate 
initially. For the disk, the radial profile of the radial ve¬ 
locity dispersion is chosen such that the Toomre Q sta¬ 
bility parameter of the 3D disk is similar to that of the 
2D counterpart considered in Paper I; we set Q « 2 over 
5 kpc < R . < 15 kpc, which is large enough to suppress a 
spontaneous growth of spiral structures when evolved in 
isolation. 

We model the perturbing companion using a dark mat¬ 
ter halo alone, which follows the truncated Hernquist 
profile 


P P tb{r) 


M p tb r pt b 
27r r(r + r pt b) 3 ’ 


for r < 100 kpc, 


(4) 


with the scale radius of r p tb = 5 kpc, and p p tb = 0 for 
r > 100 kpc. We vary the companion mass M pt b to 
study the dependence of tidal features on the strength of 
tidal force (see Section 2.2). Compared with the point- 
mass model considered in Paper I, the current extended- 
halo model realizes a more realistic, smooth variation 
of the perturbing gravitational potential especially when 
the companion is close to the pericenter. 

Most of the simulations presented in this paper use the 
galaxy model constructed by distributing N h = 4 x 10 5 , 
Nb = lx 10 4 , and N d = 1 x 10 5 particles for the halo, 
bulge, and disk of the primary galaxy, respectively, and 
TVptb = 1 x 10 5 particles for the companion. We evolve 
our 3D galaxy model in isolation for 2 Gyr, finding that 
the disk remains almost axisymmetic, with very weak (~ 
1% density variations) nonaxisymmetric features forming 
at its outskirts. This suggests that our initial disk is 
globally stable, due to a strong bulge as well as quite 
large value Q (~ 2), in the absence of tidal perturbations. 
To check the dependency of simulation outcomes on the 
particle numbers, we also run a high-resolution model 
using N h = 8 x 10 6 , N b = 2 x 10 5 , N d = 2 x 10 6 , and 
jVptb = 2 x 10 6 particles. We find that resolution does 
not make significant differences in the properties of tidal 
features: the spiral arms in the low-resolution model are 
~ 10% weaker and have a ~ 10% smaller pitch angle 
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Table 1 

Model parameters and simulation results 


Model 3, 

(1) 

-^ptb/ Mg 
(2) 

Rperi (kpc) 

(3) 

P 

(4) 

S 

(5) 

^tail 

(6) 

tan ^taii 

(7) 

^tail/^20 

(8) 

max 

(9) 

TA1 

0.44 

25 

0.410 

0.248 

0.12 

0.69 

30.9 

0.173 

TA2 

0.44 

35 

0.180 

0.166 

0.16 

0.61 

23.1 

0.108 

TA2H 

0.44 

35 

0.175 

0.163 

0.15 

0.61 

25.5 

0.113 

TA3 

0.44 

45 

0.096 

0.121 

0.21 

0.55 

13.6 

0.066 

TB1 

0.22 

25 

0.202 

0.134 

0.19 

0.45 

14.2 

0.111 

TB2 

0.22 

35 

0.089 

0.091 

0.25 

0.45 

10.2 

0.061 

TB3 

0.22 

45 

0.046 

0.065 

0.35 

0.39 

7.6 

0.039 


a Model TA2H is identical to model TA2 except that the former employs 20 times more particles than 
the latter. 


than in the high-resolution counterpart. 


2.2. Model Parameters 


Initially, we place the two galaxies at large separation, 
and make them move on mutual parabolic orbits with 
a pericenter distance R pcr i that is calculated under the 
assumption that the total mass of each galaxy is con¬ 
centrated at its own center of mass. The galaxy orbits 
are prograde and located at the same plane as the disk 
of the primary galaxy. To study tidal interactions with 
differing strength, we consider six models that differ in 
T/ptb and R per [. 

The responses of the stellar disk to the tidal field of 
the companion can be parameterized by either 


P = 


Tfptb 


Ra 


r p tb 


R 


peri 


(5) 


or 



R n 


1'pt.b T Rp 



( 6 ) 


where R g , taken equal to 7 R^ in this work, is the char¬ 
acteristic size of the primary galaxy (e.g., Springel & 
White 1999), M g is the galaxy mass within R g , AT is 
time elapsed for the companion to orbit one radian near 
the pericenter, and T is time taken by the disk particles 
at R = R g to rotate one radian about the galaxy center. 
The P parameter in Equation (5) directly measures the 
tidal force exerted by the companion relative to the grav¬ 
itational force of the primary at R g (e.g., Byrd & Howard 
1992). On the other hand, the S parameter in Equation 
(6) takes allowance for the duration of tidal interaction 
and thus corresponds to the ratio of the linear momen¬ 
tum imparted to a disk particle by the companion with 
respect to its original linear momentum of galaxy rota¬ 
tion (e.g., Elmegreen et al. 1991). Note that the current 
definition of S differs slightly from that given in Paper I, 
since the companion now has an extended density distri¬ 
bution rather than being treated as a point mass. When 
the tidal perturbations are applied impulsively over the 
time scale of AT, S' is a good measure of tidal strength 
(e.g., Section 8.2 of Binney & Tremaine 2008). One can 
thus expect that S is a better measure of tidal forcing 
for tidal tails formed in the outer disk where the rotation 
period is longer than AT. 

Table 1 lists the parameters of each model and some 
simulation outcomes. Column (1) labels each run. We 
choose model TA2 with M pt h/M g = 0.44 and R per i = 
35 kpc as our fiducial model, which is the 3D counter¬ 
part of model A2 of Paper I. Model TA2H is identical to 


model TA2 except that the former employs 20 times more 
particles than the latter. Column (2) gives the compan¬ 
ion mass relative to the primary galaxy. Column (3) lists 
the pericenter distance of two galaxies under the Keple- 
rian orbits. Columns (4) and (5) give the dimensionless 
interaction strengths P and S', respectively. Columns 
(6), (7), and (8) list the time f ta ii when the tidal tail de¬ 
velops strongest, and its pitch angle tani ta ii and surface 
density E ta ii/E 20 at t = t ta u, respectively. Here, E 20 is 
the surface density of the initial disk at R = 20 kpc. Fi¬ 
nally, Column (9) gives the peak strength of the spiral 
arms that are induced. 

We take 10 10 M 0 , 1 kpc, and 1 km s -1 as the units of 
mass, length, and velocity, respectively. In these units, 
the gravitational constant is G = 4.289x 10 4 , and the cor¬ 
responding unit of time is t 0 = 0.98 Gyr. The simulation 
time is set such that t = 0 corresponds to the closest ap¬ 
proach of the galaxies under the prescribed parabolic or¬ 
bits, but the tidal deformation of the galaxies makes the 
actual time of the closest approach delayed to t = 0.02- 
0.04, with a smaller value corresponding to weaker tidal 
forcing. All the simulations start from t = —1.5, when 
the two galaxies are separated more than 280 kpc. For 
model TA2, the initial positions and velocities of the 
galaxies are rq = {R\,6\,z\) = (80.3,—2.4, 0) and rq = 
(—39.0,0.2,0) for the primary, and r 2 = (206.3,0.7,0) 
and r 2 = (—100.2,0.2,0) for the companion. 

The simulations are performed using the GADGET 
code that employs the Barnes-Hut hierarchical tree algo¬ 
rithm (Barnes & Hut 1986) to solve the Poisson equation 
(Springel et al. 2001). We adopt the force error tolerances 
of a = 0.02 and 9 = 0.8. The gravitational softening pa¬ 
rameter in the low-resolution models is set to 0.1 kpc for 
the disk, the halo, and the companion, and 0.07 kpc for 
the bulge, corresponding to 2.8 times the mean particle 
distance in each system (e.g. Springel et al. 2001). 

3. PROPERTIES OF TIDAL FEATURES 

Tidal perturbations disturb not only the otherwise ax- 
isymmetric disk, producing a tidal tail and a bridge in 
the outer parts (R > 20 kpc) and spiral waves in the in¬ 
ner parts (4 < R < 15 kpc), but also the primary halo 
and companion galaxy. In this section, we explore the 
correlations between the physical properties of the tidal 
features and the interaction strength parameters P and 
S. The deformation of the halo and companion will be 
presented in Section 5. 

3.1. Tidal Tail and Bridge 

Figure 2 shows morphological changes of the stellar 
disk of the primary galaxy in our fiducial model TA2 
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Figure 2. Distributions of the disk particles of the primary galaxy in model TA2 projected on the X—Y plane. In each panel, the position 
of the companion is marked by a grey circle, with its size indicating the half-mass radius. The solid lines draw the actual trajectories of 
the primary and companion, while the dashed lines in the top-left panel give their Keplerian orbits that assume no change in the mass 
distributions of the two galaxies. 
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projected on to the orbital plane. The grey circle in each 
frame indicates the center of mass of the companion at 
a given epoch, with its size corresponding to the region 
encompassing a half of the companion mass. The solid 
curves represent the actual trajectories of the two galax¬ 
ies. The dashed lines in the t = —0.5 frame draw the 
parabolic orbits they would follow if their extended mass 
distributions were undisturbed over the course of the en¬ 
counter. The difference between the solid and dashed 
orbits are due to the angular momentum loss caused by 
dynamical friction, as will be discussed in Section 5 in 
detail. 

At early time when the galaxies are separated widely 
(t < —0.1), tidal perturbations are so weak that the 
orbits are parabolic and the disk of the primary re¬ 
mains almost axisymmetric. As the perturbing galaxy 
approaches the pericenter, the tidal force begins to ex¬ 
cite the epicycle motions of stars in the outer disk, first 
shaping them into a bridge at the near side to the com¬ 
panion (t ~ 0 — 0.2) and then a tail at the opposite side 
(t ~ 0.1 — 0.3). In the bridge, the epicycle phases of the 
constituting particles are forced to be aligned to the com¬ 
panion in such a manner that the radial velocities and 
the gradient of the azimuthal velocities are always max¬ 
imized along the line connecting the galaxies (Paper I). 
On the other hand, the tidal tail develops as strongly 
perturbed particles at the near side catch up and overlap 
with mildly perturbed particles at the far side (Paper I; 
Pfleiderer 1963; Toomre & Toomre 1972). The tail is 
approximately logarithmic in shape, indicating that its 
pitch angle is independent of R. 

The tails and bridges produced in our models are tran¬ 
sient and decay after t ~ 0.2 — 0.4. As the companion re¬ 
cedes, the tidal force becomes weaker and is thus unable 
to keep the coherency of the forced epicycle phases of the 
bridge particles. Composed of particles gathered from a 
wide range of radii in the unperturbed disk, the tails have 
a large velocity dispersion, tending to disperse with time. 
They disappear almost completely at i ~ 1.0. A small 
fraction (~ 2%) of the disk particles that achieve veloc¬ 
ities larger than the escape velocity become unbound, 
either captured by the companion or escaping from the 
system. With no further perturbations from the compan¬ 
ion at t > 1.0, on the other hand, the bound particles 
once in the bridge or tail spread out, following galaxy 
rotation with large eccentricities. The disk becomes in¬ 
creasingly featureless over time. 

Paper I found that the formation time, shape, and 
strength of the tails in 2D models are well correlated 
with S. To study how the correlations change in the 
current self-consistent 3D models, we define the tail for¬ 
mation epoch, ttaih as the time when it becomes dens¬ 
est at R = 20 kpc, and measure its pitch angle 'bail 
and surface density E tai] at that time. Figure 3 plot 
these values as filled symbols against P and S ; these are 
also tabulated in Columns (6)-(8) of Table 1. With 
weak tidal interaction, model TB3 does not produce a 
readily identifiable tail, suggesting that the formation 
of tidal tail requires P > 0.05 or S > 0.07. This 
threshold is probably a lower limit since all of our mod¬ 
els consider in-plane, prograde encounters, the most fa¬ 
vorable condition for the tail development. A stronger 
tail tends to form earlier and has a larger pitch an¬ 
gle. The dashed lines are the best power-law fits to the 


3D results: f t aii = 0.07P 0 49 , tanz ta ii = 0.1P 0 24 , and 
Staii/S 2 o = 60P 0 65 in terms of P; i tBil = 0.04S- 079 , 
tanitaii = 1.25S' 0 43 , and E t aii/E 2 o = 1 ^S 1A in terms of 
S. As expected, the tail properties are better correlated 
with S than P: the linear fitting coefficients are ~ 0.82 
and ~ 0.93 against P and S, respectively. For compari¬ 
son, we also plot the 2D results taken from Paper I as as 
open circles. Overall, the 2D and 3D results agree very 
well, although the 3D tails are slightly weaker, which is 
probably caused by weaker disk gravity due to finite disk 
thickness. Figure 3 shows that the tail in model TA2H 
with high resolution has almost the same properties as 
in model TA2, suggesting that our results do not depend 
sensitively on resolution. 

Another way to quantify the tail strength is to use 
the tidal response introduced by Springel & White 
(1999) who ran a number of iV-body simulations for 
galaxy mergers with differing halo mass and interaction 
strength. They defined 7~ e g as the maximum fraction 
of the disk particles in each model that reach distances 
beyond lOP^ (= 34 kpc in our models) from the disk 
center in the course of interaction. Figure 4(a) plots %g 
from our tail-forming models (filled circles) as well as the 
models (open symbols) from Springel & White (1999) as 
a function of S: regular triangles are for models with 
Rperi/Rd > 7; inverted triangles for 2.5 < R pe ri/Rd < 7; 
squares for R per i/Rd < 2.5. All of our models have 
Pperi/Pd > 7. Despite differences in the galaxy mod¬ 
els such as potential depth and rotation curve, etc., our 
results are overall in good agreement with the extrapo¬ 
lation of Springel & White (1999) for S < 1. Note that 
in very strong encounter models with S' > 1 in Springel 
& White (1999) represented by open squares, the com¬ 
panion directly passes through the inner parts of the pri¬ 
mary disk, for which the S parameter relying on a dis¬ 
tant tide approximation is not useful in characterizing 
the interaction strength. Except for these very strong- 
encounter models, there is a good correlation between 
T e g and S with some scatters. The best fit for models 
with Rperi/Rd > 2.5 is T e g = 0.16S 1 - 1 , which is shown 
as the solid line in Figure 4(a). 

Since Tcfi includes the particles not only in a tail but 
also escaped or captured to the companion, to mea¬ 
sure the tail strength we similarly define 7taii as the 
fraction of the disk particles consisting only of a tail 
at R > 10 Rd- Figure 4(b) plots "7t a ii for our models 
as a function of S. The solid line draws the best fit 
7taii = 0.315 121 , which has a tighter correlation than 
the T e g-S relation. The 7t a ii _ 5' relationship can be ex¬ 
plained as follows. Let us assume that the tidal perturba¬ 
tions are applied impulsively near the pericenter during 
the time interval AT. The radial velocity increment at 
R g is then given by Svr = 2GM p R g AT /(i? per i + r p ) 3 . 
This enhances the epicycle amplitudes of the disk par¬ 
ticles to SR ~ Si’u/Kg. With the epicycle frequency 
K g = \/2n g = (2GMg/i? 2 ) 1 / 2 at R g for fiat rotation, 
it follows that SR/R g = \f2S. Since the perturbed mass 
that goes into the tail is proportional to T, 0 SR 2 , one ob¬ 
tains Ttaii oc S 2 , similar to our numerical results. 

3.2. Spiral Arms 

We have seen in the previous section that the tidal 
force from the companion provides strong perturbations 
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Figure 3. Dependence of the formation epoch t ta ii, the pith angle Fa.il ? and the surface density S ta ;i of tidal tails on the tidal strength 
parameters P (left) and S (right). Filled symbols are the current 3D results, while the 2D razor-thin results are compared as open circles. 
The dashed lines are the best fits: t ta ii = 0.07P - °- 49 , tan F.a.il = 0.1P 0 ' 24 , and Eta.il/^20 = 60P 0 ' 65 in terms of P; f ta ii = 0.04S - °' 79 , 
tani ta ;i = 1.25S 0,43 , and S ta ;i/S 2 o = 1445 11 in terms of S. 


to epicycle orbits of particles in the outer disk, resulting 
in a tidal tail and a bridge. Tidal disturbances in the 
inner part of the disk are not as strong as in the outer 
parts, but can nevertheless induce spiral arms there (e.g., 
Toomre 1969; Donner et al. 1991). Figure 5 plots close- 
up snapshots of the stellar surface density in logarithmic 
scale in model TA2H projected onto the orbital plane, 
with the center of mass of the disk shifted to the origin 
(SX = SY = 0). At t = 0.1, the outer edges of the arms 
are quite sharp, which is a common feature of tidally- 
induced arms (e.g., Struck 1990; Elmegreen et al. 1991; 
Elmegreen et al. 2011; Dobbs et al. 2010). 

To delineates the spiral structure, Figure 6 plots the 
perturbed surface density, <S£ = £ — E 0 , relative to the 
initial surface density £o in the pin R plane. Before the 
companion passes by the pericenter (t < —0.1), the disk 
appears almost axisymmetric, without forming a bar or 
spiral structures, indicating again that the disk is glob¬ 
ally stable in isolation. Tidal disturbances excite epicycle 
motions of the particles in the inner disk, which are co¬ 
herently organized to build up a well-defined two-armed 
global spiral pattern. With Q ~ 2, the effect of swing 
amplification on the growth of the spiral arms is not sig¬ 
nificant (Paper I). It is apparent that the arms are ap¬ 
proximately logarithmic. In model TA2H, they become 
strongest at t = 0.2 —0.4 and decay as they keep winding 
out over time. In these single encounter experiments, no 


spiral arms are apparent after t = 1.5 (see also Fig. 5). 

3.2.1. Arm Strength 

One of the most important parameters that control the 
responses of gas flows across a stellar spiral arm is the 
arm strength (e.g., Kim & Kim 2014; Kim et al. 2014). 
It has often been customary to quantify the arm strength 
using the dimensionless parameter 

_ 27rG£ m= 2 , , 

Rn 2 ’ 1 ’ 

where £ m= 2 denotes the amplitude of the m = 2 Fourier 
mode of £. Note that J- effectively measures the maxi¬ 
mum gravitational force due to the spiral arms as a frac¬ 
tion of the axisymmetric gravitational force in the un¬ 
perturbed disk (Roberts 1969; Shu et al. 1973; Kim & 
Ostriker 2002, 2006; Shetty & Ostriker 2006). 6 

Figure 7(a) plots the radial variation of T averaged 
over the time interval At = 0.45 centered at the epoch 

6 For m- armed tightly-wound spirals with pitch angle i(<g 1) and 
radial wavenumber k,R = m/(R tan i), the corresponding gravita¬ 
tional potential at z = 0 of a razor-thin disk is = <E> m exp 
ikftR) with <l> m = —2nGT l m/kji (e.g., Eq. (6.31) of Binney 
& Tremaine 2008). Therefore, T = \d&m/dR\ max /(RQ 2 ) = 
m|<l m |/(tanlR 2 fl 2 ) = 27rGE m /(Ml 2 ). 
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Figure 4. Dependence on S of (a) 7^ff, the maximum fraction 
of the disk particles at R > 10 Rd, and (b) 7tail> the fraction of 
the disk particles in the tail at R > 10 Rd- Filled circles are from 
our 3D simulations. The results of Springel Sc White (1999) for 
models with R pe ri/Rd > 7 (regular triangles), 2.5 < R pe ri/Rd < 7 
(inverted triangles), and R pe ri/Rd < 2.5 (squares) are compared. 
The solid lines are the best power-law fits: 7^ff = 0.165 11 and 
Ttail = 0.31S 2 - 1 . 

when the arms attain the maximum strength in each 
model, while Figure 7(b) gives the temporal change of 
T averaged over a range of radii where the arms are 
strong. Definitely, a weaker tidal interaction produces 
weaker arms that grow more slowly at larger galactocen- 
tric radii: they are maximized at R ~ 5 — 10 kpc for 
the intermediate encounter models TAl, TA2, and TB1 
with P > 0.17 (or S > 0.13), and at R ~ 8-15 kpc for 
the weak encounter models TA3, TB2, and TB3 with 
P < 0.1 (or S < 0.12). The difference between the re¬ 
sults of models TA2 and TA2H is insignificant, indicating 
again that resolution does not affect the simulation out¬ 
comes much. For comparison, Figure 7 also plots the 
results of the 2D model A2 taken from Paper I as dashed 
lines. Note that the arms in model A2 have T about 
70% stronger than in the 3D counterpart, which is pri¬ 
marily due to the overestimated self-gravity in the 2D, 
razor-thin geometry. 

Figure 8 plots the peak arm strength .F max as func¬ 
tions of P and S. The solid lines represent the best 
power-law fits P max = 0.39P 07 and P max = 0.82S 1 ' 1 to 
the 3D results, with the linear correlation coefficients of 
0.99 and 0.95 relative to P and S, respectively. These 
can be compared to the 2D results of P max = 0.79S' 0 83 
from Paper I shown as dotted lines. Observations usu¬ 
ally measure the arm strength in terms of the relative 
Fourier amplitude TZ = £ m= 2 /£ m= o (e.g., Rix & Rieke 
1993; Patsis, Heraudeau, & Grosbpl 2001; Kendall et al. 
2015). In order to facilitate comparisons with observa¬ 
tions, Figure 9 presents the dependence on P of TZ cal¬ 
culated from our 3D results. The filled circles plot the 
averages over 0 < t < 1 and 5 kpc < R < 15 kpc, with 
the errorbars corresponding to the standard deviations. 
The solid line is the best fit: TZ = 1.22 + 0.70 log P. 


The spiral arms achieve the full strength at t ss 0.3-0.6 
with smaller values corresponding to stronger tidal inter¬ 
action, after which they begin to decay. This is mainly 
caused by the increase in the velocity dispersions of the 
disk particles. First of all, the growth of spiral arms re¬ 
quires the particles at different radii to gather into the 
arms, which increases the velocity dispersions (Paper I). 
Also, heating of the disk particles due to gravitational 
scattering off the arms becomes efficient once the arms 
are strong enough, inhibiting further growth of the arms 
(e.g., Sellwood & Carlberg 1984; Binney 2001). The in¬ 
creased velocity dispersions make the epicycle orbits in 
the arms kinematically less coherent. Figure 7(b) shows 
that the amplitudes of the arms in the current 3D models 
decrease with time almost exponentially, with a charac¬ 
teristic timescale of ~ 0.5 Gyr, about twice faster than in 
the 2D models of Paper I. The faster decay of the arms in 
the 3D models is presumably due to weaker disk gravity 
in a stratified disk, which makes the arms wind out at a 
faster rate as well. 

3.2.2. Pitch Angle 

Figure 6 shows that the tidally-induced spiral arms are 
approximately logarithmic in shape over a wide range of 
radii. This appears consistent with the observational re¬ 
sults that grand design spirals are close to be logarithmic 
in optical and near-infrared images (e.g., Grosbpl & Pat- 
sis 1998; Seigar et al. 2006; Davis et al. 2012; Martinez- 
Garcfa 2012). As in Paper I, we calculate the arm pitch 
angle at a given time by first finding the Fourier coeffi¬ 
cients in <j> and In R defined as 

a (p) = j-j~ 5Z CX P[*( 2< A? + P lnR j )]> (8) 

p j —i 

where N p is the number of particles located at R = 5- 
10 kpc for the intermediate tidal-forcing models with 
P > 0.17 (or S > 0.13), and at R = 8-15 kpc for the 
weak forcing models with P < 0.1 (or S < 0.12), ( Rj,(pj ) 
are the coordinates of the j- th particle, and p is a real 
number corresponding to the slope of a two-armed spiral 
in the (Ini?, 0) plane. We then calculate p max that max¬ 
imizes |A(p)| and calculate the pitch angle of the arms 
through tan* = 2/p max at each time (e.g., Sellwood & 
Carlberg 1984; Sellwood & Athanassoula 1986). 

Figure 10 plots the temporal changes of the arm pitch 
angles for all 3D models. The 2D result of model A2 
from Paper I is also compared. When the arms stand 
out initially, they have moderate pitch angles amount¬ 
ing to i ~ 15° — 20°. After attaining the peak ampli¬ 
tude, the arms in the 3D models wind out over time as 
tan* oc t~ 0 - 75 , which is steeper than tan* oc t~ 06 in 
the 2D models. The faster winding rates in 3D models 
are again because of reduced self-gravity in vertically- 
extended disks. Note that the pitch angle of purely 
kinematic density waves with no self-gravity decays as 
tan* oc t^ 1 (e.g., Section 6.2 of Binney & Tremaine 
2008). 

3.2.3. Pattern Speed 

The pattern speed of spiral arms is also an important 
parameter, although it is not well constrained observa- 
tionally. To calculate the arm pattern speed induced in 
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Figure 5. Zoom-in snapshots of the disk surface density in model TA2H at the times shown in Figure 2. The coordinates are shifted so 
as to make the center of mass of the disk located at the origin (<5X = SY = 0). The colorbar labels log(E/10 4 Mqpc - 2 ). The spiral arms 
become strongest at 0.3 — 0.4 and decay afterwards. 
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Figure 6. Distributions of the perturbed surface density S S/So of model TA2H in the In R plane. When t < 0, the disk is almost 
featureless. When 0 < t < 0.2, S S is dominated by a bridge and a tail at R > 15 kpc, which become weak after t ~ 0.3. The spiral arms at 
R < 15 kpc are approximately logarithmic and wind progressively with time. The greyscale bar labels log<5E/Eo. 
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t 

Figure 7. (a) Radial dependence of the arm strength T averaged 

over the time interval At = 0.45 centered at the time of the peak 
strength, (b) Temporal dependence of T averaged over 5 kpc < 

R < 10 kpc for models TA1, TA2, and TB1, and over 8 kpc < R < 
15 kpc for models TA3, TB2, and TB3. After the peak, T decays 
as ~ exp(—1/0.5 Gyr). In both panels, the results of the 2D model 
A2 from Paper I are compared as dashed lines. 

our models, we use the normalized cross-correlation of 
the perturbed surface densities at two different epoches 
separated by At as 

1 f 2n 

C{R,<p,t) = / 6T,(R,<j),t)ST,(R,<j>+ip,t+At)d(t). 

So [R)~ Jo 

(9) 

By taking a sufficiently small value of At = 0.1, we 
find </ 2 max that maximizes C(R , ip,t) at given radius and 
time. The instantaneous pattern speed is then given by 
t ) = (jCmax/At. 

Figure 11 plots as contours the amplitudes of the 
C(R,<p,t) in the R-(ip/At) plane for some selected 
epoches of models TA2 (left) and TB2 (right). The an¬ 
gular frequencies Q and fi± n/2 from the initial rotation 
curve are overplotted as a solid line and two dotted lines, 
respectively. The locus of </3 max hr each panel draws the 
arm pattern speed as a function of R. Shortly after the 
pericenter passage of the companion (t ~ 0.1), the per¬ 
turbed surface density is dominated by the tidal bridge 
rather than the arms. The bridge has an almost fixed 
pattern speed at S2 p ~ 10 km s -1 kpc -1 , which corre¬ 
sponds to the angular speed of the companion near the 
pericenter. At this time, the particles inside the bridge 
are tidally locked to the companion. As the compan¬ 
ion moves away from the primary, the bridge becomes 
weaker and C(R,ip,t) becomes progressively dominated 
by the spiral arms. At f ~ 0.3 when the arms have sub¬ 
stantial amplitudes, 52 p is larger than f1 — n/ 2, but only 
slightly. Since Q p is not constant over radius, the arms 
produced in our models are not exactly a pattern in a 
strict sense. They are rather similar to gravity-modified 
kinematic density waves. Although self-gravity tends to 
increase 52 p , its effect is not significant in our 3D models. 
At late time (t > 0.6), Q p converges to the 52 — k/2 curve. 

To summarize this section, spiral arms induced by a 


tidal interaction in our models are not quasi-stationary 
density waves envisaged by Lin & Shu (1964). They are 
rather kinematic density waves slightly modified by grav¬ 
ity. Due to reduced self-gravity, spiral arms in the 3D 
models are weaker and decay faster than those in the 
razor-thin counterparts. The arms in vertically-stratified 
disks have a smaller pitch angle and wind up more rapidly 
than those in the razor-thin, 2D disks. 


4. DENSITY STRUCTURE OF SPIRAL ARMS 

While there are numerous studies on the generation 
of stellar spiral arms, they concentrate mostly on arm 
morphologies and longevity without focusing on the arm 
density structures (e.g., Hernquist 1990b; Barnes 1992; 
Mihos & Hernquist 1994; Naab & Burkert 2003). With¬ 
out much information on the vertical structure of stellar 
spiral arms, most previous works that studied galactic 
spiral shocks across the arms employed the arm poten¬ 
tials that are independent of, or varying weakly with, 
the vertical height (e.g., Kim & Ostriker 2006; Kim et 
al. 2006, 2010). Cox & Gomez (2002) suggested an an¬ 
alytic expression for 3D spiral density perturbations as 
described below. In this section, we analyze the den¬ 
sity structures of spiral arms produced in our models 
and compare them with the analytic suggestion of Cox 
& Gomez (2002). 

Starting from a physically-motivated trial function for 
the spiral gravitational potential that drops off exponen¬ 
tially in the radial direction, Cox & Gomez (2002) ob¬ 
tained a simple formula for the m = 2 spiral density 
perturbation 


Pi(R><t>,z) = Poexp - 


R R, 


with 


7 = 2 


R, 


(f) 0p(i^o) 


n—1 


n cos(nj), 

MR/RoY 


tan i 


( 10 ) 

( 11 ) 


where po is the normalization coefficient, i?o is the fidu¬ 
cial radius, R s is the radial scale length of the perturbed 
density, h is the vertical scale height, a n is the Fourier 
amplitude of an azimuthal mode n , and <fi p (Ro) denotes 
the azimuthal phase of the arm at R = Ro. The sum¬ 
mation over n is to allow for linear superpositions of 
various modes along the azimuthal direction for two¬ 
armed, logarithmic spirals. As examples, Cox & Gomez 
(2002) considered two different cases: sinusoidal arms 
having N = 1 with oi = 1 (e.g. Kim & Ostriker 2002; 
Wada & Koda 2004; Shetty & Ostriker 2006) and concen¬ 
trated arms having N = 3 with ai = 8/(37 t), <22 = 1/2, 
and ai = 8/(157r). The concentrated arms have flatter 
interarm regions than the sinusoidal arms. The four¬ 
armed concentrated spirals were used by Dobbs & Bon- 
ncll (2006) and Dobbs et al. (2006) to study gas dynamics 
driven by spiral arms. The thick spiral potential con¬ 
sidered by Patsis, & Grosbpl (1996) corresponds to the 
sinusoidal arms. 

Figure 12 plots the vertical density profiles in the arm 
regions at R = 6 kpc and <j> = 90° of model TA2H at t = 
0.38 when the arm strength is nearly peaked (left), and at 
t = 0.85 when the arms are in the decaying phase (right). 
The black solid, dotted, and dot-dashed lines draw the 
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Figure 8. Dependence of the peak arm strength J ma x on (a) P and (b) S. A stronger tidal forcing results in stronger arms. The solid 
lines are the best fits, J max = 0.39P 0 7 and P m ax = 0.82S' 11 , to the 3D results, while the dotted lines draw the 2D results adopted from 
Paper I. 



Figure 9. Dependence of the relative amplitude 1Z = 
S m = 2 /S m= o of the arms on the P parameter. The solid circles 
mark the averaged over 0 < t < 1 and 5 kpc < R < 15 kpc, with 
the errorbars indicating the standard deviations. The solid line is 
the best fit: 1Z = 1.22 + 0.70 log P. 

total density p(R, <j>, z), the azimuthally-averaged density 
p(R,z) = f p(R, 4>, z)d(/)/2ir, and the perturbed density 
Pi = p — p 1 respectively. The three density distributions 
are well described by sech 2 (z/h), shown as thin red lines, 
with h = 0.33, h = 0.36, hi = 0.30 kpc at t = 0.38 for 
the total, azimuthally-averaged,_perturbed densities, re¬ 
spectively, and with h = 0.37, h = 0.36, hi = 0.41 kpc 
at t = 0.85. The temporal increase of the scale heights 
between t = 0.38 and 0.85 is due to disk heating. But, 
these agree within 10%, validating the vertical depen¬ 
dence of pi in Equation (10), with a scale height similar 
to that of the background averaged disk. 



Figure 10. Temporal variations of the pitch angle i of the spiral 
arms located at 5 kpc < R < 10 kpc for the intermediate forcing 
models TA1, TA2, and TB1, and at 8 kpc < R < 15 kpc for the 
weak forcing models TA3, TB2, and TB3. The result of the 2D, 
razor-thin model A2 is compared as a thick grey line. The arm 
pitch angle in the 3D models decays as t ~ 0,75 after the peak, which 
is steeper than /. IMi in the 2D model. 

Figure 13 plots as solid lines pi as a function of R along 
the cuts with </> = ±90° for model TA2H in the mid¬ 
plane. The top and second rows are at t = 0.38, while 
the third and bottom rows are at t = 0.85. The thin solid 
lines give the simulation results, while the thick grey lines 
draw the best fits using Equation (10) based on the sinu¬ 
soidal arms (left panels) and the concentrated arms (right 
panels). The fact that Equation (10) describes the arm 
positions represented by the over-densities fairly well im- 
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Figure 11. Contours of the cross-correlation C(R,ip,t.) of the perturbed surface density in the radius-frequency plane for models TA2 
(left) and TB2 (right) at t = 0.1, 0.2, 0.3, 0.6, and 0.9. The solid and dashed lines draw 12 and 12 ± re/2 curves. At t < 0.2, the cross¬ 
correlation is dominated by the tidal bridge and tail at R > 17 kpc, while the spiral arms dominate at 5 kpc < R < 15 kpc for t > 0.3. 
Note that 12 p of the arms is very close to the 12 — re/2 curve. 


plies that the spiral arms are closely logarithmic in shape. 
Overall, the concentrated arms fit the radial arm struc¬ 
ture better at t = 0.38 when the arms are strong with 
T > 10%, whereas the sinusoidal arms provide a better 
fit at t = 0.85 when the arms are weak with T < 10%. 
We note however that these fits are reasonably good only 
at R > 4 kpc and fail at smaller R. This is expected from 
the fact that the tidally-driven arms do not decline mono- 
tonically with R but are peaked at 5 kpc < R < 10 kpc 
(see Fig. 7), while Equation (10) requires arms to decay 
exponentially in the radial direction over all radii. 

5. ORBITAL DECAY 


As Figure 2 shows, the actual trajectories of the pri¬ 
mary and companion galaxies during the tidal encounter 
deviate from those under the assumption of the rigid ha¬ 
los. The discrepancy is due to dynamical friction occur¬ 
ring when two galaxies overlap at least partly. As the 
companion enters the primary halo, the former creates 
a density wake in the latter that in turn exerts gravita¬ 
tional drag force on the former. The drag force removes 
the orbital angular momentum from the companion, de¬ 
creasing the pericenter distance. The orbital angular 
momentum extracted from the companion is transferred 
to the spin angular momentum of the primary halo. A 
similar process also occurs when the primary enters the 
companion halo, which results in the angular momentum 
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Figure 12. Vertical distributions of the stellar surface density in the arm regions of model TA2H at R = 6 kpc for (a) t = 0.38 when 
the arms are strong, and (b) t = 0.85 when the arms become weak. The black solid, dotted, dot-dashed lines draw the total density p, the 
azimuthally-averaged density p, and the perturbed density p\=p p. respectively. The thin red lines are fits based on the sech 2 {z/h) 
function, with the scale heights h = 0.33, h = 0.36, hi = 0.30 kpc at t = 0.38 and h = 0.37, h = 0.36, h\ = 0.41 kpc at t = 0.85, for the 
total, azimuthally-averaged, perturbed densities, respectively. 


transfer from the orbit of the primary to the spin of the 
companion. 

To quantify the amount of the angular momentum 
transfer, Figure 14 plots the temporal variations of the z- 
component of the angular momenta of all the components 
as functions of time. In each panel, the solid line draws 
the net (orbital plus spin) angular momentum, while the 
dotted and dashed lines give the orbital and spin angular 
momenta defined, respectively, by 

Lorb = MX. CM X V CM: (12) 

and 

Lspin = ^ ' U7.j(x ? ; Xqm) X (v.j V(7 m)j (13) 


where M = ^rm is the total mass, Xcm = JJmjXj/M 
is the position vector of the center of mass, and Vcm = 
JO m.iVi/M is the mean velocity of each component, mea¬ 
sured in the center-of-mass frame of the whole system. 
The dot-dashed line in Figure 14(a) plots the total an¬ 
gular momentum of the whole system, showing that it 
is conserved within 3% throughout the entire evolution. 
Initially, all the components except the disk do not spin, 
and their angular momenta are dominated by the or¬ 
bital motions. As the primary and companion approach 
the pericenter, they start to experience dynamical fric¬ 
tion, and the primary halo and the companion begin to 
continuously acquire the spin angular momentum at the 
expense of the orbital angular momentum. 

In order to check if the decrease in the orbital angular 
momentum of the companion is really caused by dynam¬ 
ical friction due to the primary halo, we use the Chan¬ 
drasekhar (1943) formula for the drag force 


Fdf = - 


47 TG 2 M^ th p h In A 


ptb 


erf(A') — ^=e x 

\/7T 


p{tJ,4) 


where ph is the average density of the halo, v p tb is the 
relative speed of the companion, A" = u p tb/(v / 2c r ), with 
a denoting the mean velocity dispersion of halo particles, 
and In A = ln(& max /6 m i n ) is the Coulomb logarithm in¬ 
troduced to avoid a singularity in the force evaluation, 
with & max and b m - ln representing the maximum and mini¬ 
mum impact parameters of the background particles, re¬ 
spectively. For orbital decay of a satellite galaxy in a 
dark matter halo, Equation (14) matches the numerical 
results provided b max is taken to be of order of the or¬ 
bital radius of the halo, and 6 m i n ss max(GM pt b/r , p tbJ r h ) 
with rh being the half-mass radius (e.g., Lin & Tremaine 
1983; Weinberg 1986). To apply Equation (14) to model 
TA2H at t = 0, we calculate the mean density of the 
halo as ph = f p h (\r-r h \)p pth (\r-r pth \)d 3 r/ f p p tb (|r- 
r p tb|)d 3 F = 3.0 x 10 5 M 0 kpc -3 , where and r pt b de¬ 
note the position vectors of the halo and the compan¬ 
ion, respectively. The velocity of the primary relative to 
the halo is u pt b = 380 km s -1 . With 6 max = R pcr ; and 
&min = (1 + V / 2)r p tb, In A = 1.3. The resulting torque on 
the companion is given by 


dL, 


'orb,ptb 

dt 


= ~ Rn 


x F 


DF, 


(15) 


which is plotted as a thick solid line in Figure 14(a), in 
fairly good agreement with the instantaneous decreasing 
rate of the orbital angular momentum of the companion. 
This validates that dynamical friction is indeed a primary 
cause of the angular momentum transfer in our numerical 
simulations. 

A close inspection of Figure 14 reveals that while the 
orbital angular momenta of the halo and companion de¬ 
crease monotonically with time, those of the disk and 
bulge stay almost constant at t < —0.1, increase slightly 
for —0.1 < t < 0.15, and subsequently decrease slowly. 
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Figure 13. Perturbed density pi (R. <f>, z) at the 2 = 0 plane as a function of R for two different azimuthal directions (o = ±90°) of model 
TA2H at t = 0.38 when the spiral arms are quite strong (top and second rows) and at t = 0.85 when the arms are weak (third and bottom 
rows). The thick grey lines represent the best fits using Equation (10) based on the sinusoidal-arm (left column) and the concentrated-arm 
(right column). 


This indicates that the disk and bulge experience positive 
torque temporarily near the pericenter, while the halo 
and companion always receive negative torque. To ana¬ 
lyze what causes this unexpected behavior of the angular 
momenta of the disk and bulge, Figure 15 plots as dashed 
and dotted lines the torques exerted on the disk/bulge 
due to particles in the companion and the halo, respec¬ 
tively. It is apparent that the torque by the companion is 
overall positive, while the halo gives a negative torque to 
the disk/bulge. The resultant net torque, shown as the 
solid line, is positive for —0.1 < t < 0.15, which boosts 
the orbital motion of the disk/bulge system. 

To understand why the companion torque is positive, 
Figure 16 plots the density distributions of the compan¬ 


ion (gray scale with black contours) and the disk and 
bulge of the primary (white contours) at t = 0 (left) and 
t = 0.25 (right) of model TA2H in the orbital plane. The 
contour levels decrease by a factor of 10°' 5 starting from 
the innermost level of 0.1 M 0 pc -3 for both companion 
and disk/bulge. Note a clear density wake in the com¬ 
panion formed at t = 0 near ( X , Y) ~ (0, 20) kpc due to 
the primary. The filled and open star symbols in cyan 
mark the centers of mass of the halo and disk/bulge of 
the primary, respectively, while the arrows in green indi¬ 
cate the velocity vectors of the disk/bulge in the center 
of mass frame of the whole system. 

The dynamical friction of the primary due to the 
companion is stronger for the (outer) halo than the 
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Figure 14. Temporal variations of the net (solid), orbital (dotted), and spin (dashed) angular momenta of (a) the companion, and the 
(b) halo, (c) disk, and (d) bulge of the primary galaxy during the tidal encounter of model TA2H. In (a), the dot-dashed line plots the total 
angular momentum of the whole system, while the thick dashed-line segment indicates the decrease of the orbital angular momentum due 
to dynamical friction expected from Chandrasekhar formula. The companion and the primary halo that do not rotate initially acquire a 
significant amount of the spin angular momentum at the end of the encounter. The orbital angular momentum of the disk/bulge experiences 
a slight boost near t = 0 and declines afterward. 


disk/bulge owing to proximity to the companion. This 
causes the halo to move slower than the rest of the pri¬ 
mary galaxy, gradually making them displaced from each 
other. At t = 0, the center of mass of the disk/bulge is 
entering the third quadrant and moving in the negative- 
Y direction, while that of the halo is still in the second 
quadrant. At this time, the position vector of the bulge, 
indicated by the blue arrow, is inclined to the line con¬ 
necting the centers of mass of the disk/bulge and the 
companion (the red dotted arrow). In such a geomet¬ 
rical configuration, the companion pulls the disk/bulge 
forward by providing a positive torque, while the halo 


tends to slow them down. It turns out that the positive 
torque is stronger than the gravitation pull by the halo at 
this time, as Figure 15 shows. The net effect is that the 
separation between the centers of mass of the disk/bulge 
and the halo increases at a faster rate. The positive 
torque becomes strongest t = 0.06 when the companion 
is closest to the disk/bulge, and then decreases as the 
companion moves away. At t = 0.25, the positive torque 
by the companion is smaller than the negative torque by 
the halo. The disk/bulge is subsequently pulled back to, 
and moves in phase with, the halo after t ~ 1. 

Does the dynamical friction between two galaxies re- 
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Figure 15. Torques on the disk/bulge exerted by its own halo 
(dotted) and the companion (dashed) in model TA2H as functions 
of time. The solid line gives the total torque. While the halo 
torque is always negative, the companion torque is positive over 
—0.1 < t < 0.6, caused by the displacement of the disk/bulge from 
the center of mass of the halo (See Figure 16). With the positive 
net torque, the disk/bulge gains angular momentum during —0.1 < 
t < 0.15. 

ally transfer the orbital angular momentum of one galaxy 
to the spin angular momentum of the other? Figure 17 
plots the temporal changes of the axially-averaged ro¬ 
tation velocities (v^) of the halo (left) and companion 
(right). The vertical arrow in each panel marks the po¬ 
sition of the companion and the primary in the left and 
right panels, respectively. Initially (t = —1.5), they are 
non-rotating, spherical, and about 280 kpc apart. As the 
center of the companion (halo) moves through the outer 
parts of the primary halo (companion), it creates a den¬ 
sity wake in the halo (companion) that pulls it backward. 
At the same time, the wake in the halo (companion) is 
pulled forward by the companion (halo), acquiring pos¬ 
itive rotational velocities as shown in Figure 17. The 
mass-weighted mean rotational velocities of the compan¬ 
ion and the halo at t = 1.0 is ( v ^) ~ 15 km s -1 and 
7 km s -1 , respectively. This suggests that tidal interac¬ 
tions can be a source of the spin angular momentum of 
a dark matter halo that might be non-rotating when it 
first formed. 

6. SUMMARY & DISCUSSION 
6.1. Summary 

Using self-consistent, 3D iV-body simulations, we have 
investigated the physical properties of tidal features in¬ 
duced in a disk galaxy interacting with its companion. 
The disk galaxy consists of a halo, a bulge, and a stellar 
disk, corresponding to Milky-Way type galaxies, all of 
which are represented by live particles. The perturbing 
companion is simplified by a live halo alone. This work 
is a straightforward extension of Paper I that considered 
a 2D, razor-thin disk residing within a fixed halo. We 
do not consider the effect of the gaseous component in 


the present work. By varying the companion mass and 
the pericenter distance, we explore the intermediate and 
weak interaction regimes with the strength parameter of 
0-04 < P < 0.41 or 0.06 < S < 0.25 (see Eqs. [5] and [6] 
for the definitions of P and S). We analyze the properties 
of tidal tails formed in outer regions, and spiral arms in 
the inner regions, and study their dependence on P and 
S. We also study the 3D density structures of spiral arms 
as well as the orbital decay of the galaxies caused by dy¬ 
namical friction occurring during a tidal encounter. The 
main results of this work can be summarized as follows. 

1. The tidal force from the companion excites epicy¬ 
cle motions in the outer regions of the primary, 
forming a bridge at the near side and an extended 
tail at the far side for interactions with P > 0.05 
or S > 0.07: tidal forcing weaker than this is un¬ 
able to produce a tidal tail. The tail formation 
time, pitch angle, and surface density measured at 
R = 20 kpc scale as t ta ii = 0.07P -0 49 , tani ta ii = 
0.1P 0 - 24 , and E ta ii/E 20 = 60P 0 65 in terms of P, 
and tta.ii = 0.04S' -0 79 , tani ta ii = 1.255 0 - 43 , and 
St a ii/S 2 o = 144S 111 in terms of S. This is in good 
agreement with the 2D results of Paper I, although 
the 2D tails are slightly stronger than the 3D coun¬ 
terparts, due to stronger gravity in the disk mid¬ 
plane. We also find that the tail response Ttaii de¬ 
fined by the fraction of the disk particles consisting 
of a tail at R > 10 Rd has an excellent correlation 
with S as 7t a ii = 0.31S 121 . 

2. The tidal force also excites two-armed spiral den¬ 
sity waves over a range of radii in the inner parts 
of the disk. We calculate the spatial and tempo¬ 
ral dependence of the arm strength J- defined by 
the ratio of the maximum gravitational force due 
to the m = 2 modes to the centrifugal force of 
disk rotation (Eq. [7]). For intermediate tidal in¬ 
teractions with P > 0.17 (or S > 0.13), the arms 
are strongest at R ~ 5-10 kpc, while weak interac¬ 
tion models with P < 0.1 (or S < 0.12) have the 
arms induced at R ~ 8-15 kpc. The spiral arms 
are stronger in models with larger P or 5, with 
the peak arm strength given by Pmax = 0.39P 0 7 
or P max = 0.825 1 ' 1 . The spiral arms are approx¬ 
imately logarithmic in shape, with a pitch angle 
i ~ 15° — 20° when the arms are strongest. Af¬ 
ter the peak, the arms decay with a timescale of 
~ 0.5 Gyr, and wind out as tan i oc f~ 0 75 . The 
derived arm pattern speed fl p is a decreasing func¬ 
tion of R and becomes close to the U — k/2 curve 
at late time. This suggests that the induced spiral 
arms are kinematic density waves modified weakly 
by self-gravity. Compared to the 2D counterparts, 
arms in 3D models are weaker, have a smaller pitch 
angle, and wind and decay more rapidly with time, 
due to weaker gravity. 

3. When compared with the analytic expression of 
Cox & Gomez (2002), the 3D density structure 
of tidally-induced, m = 2 arms can be reasonably 
well described by the concentrated arms when the 
arms are in the nonlinear regime with J- > 10%. 
On the other hand, arms in the linear regime with 
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Figure 16. Density distributions at the z = 0 plane of the companion (grey scale and black contours) and the disk/bulge of the primary 
(white contours) of model TA2H at (left) t = 0 and (right) t = 0.25. The filled and open star symbols in cyan mark the centers of mass 
of the halo and disk/bulge of the primary, respectively. The blue and green arrows in each panel are the position and velocity vectors of 
the disk/bulge in the center of mass frame of the whole system, respectively, while the red dotted arrow indicates the direction from the 
disk/bulge to the companion. The greyscale bar labels log(p/l Mq pc -3 ). The contours mark log(p/l Mq pc -3 ) = —1.0, —1.5, —2.0, • • • 
from inside to outside. 


P < 10% are better fitted by the sinusoidal model. 
The perturbed density in the arm regions follows 
the characteristic profile pi(z) oc sech 2 ( 2 /hi) along 
the vertical direction, with the scale height hi not 
much different from that of the initial disk. 

4. The halos of the primary and companion become 
partially overlap in the course of a tidal encounter. 
This creates gravitational density wakes in the ha¬ 
los, resulting in dynamical friction. The dynam¬ 
ical friction has a few notable effects. First, it 
transports the orbital angular momentum of one 
galaxy to the spin angular momentum of the halo 
of its companion, making the pericenter distance 
smaller than that under the rigid halos. Sec¬ 
ond, the initially rotation-free halos achieve the 
mass-weighted rotational velocities amounting to 
(v^) ~ 15kms _1 and 7kms~ 1 for the primary 
and companion, respectively. Third, the dynam¬ 
ical friction is stronger for the halo, causing off¬ 
sets between the centers of mass of the halo and 
disk/bulge of the primary. Thus, the disk/bulge 
can temporarily gain angular momentum by a posi¬ 
tive torque from the companion near the pericenter, 
although they are later pulled back to the center of 
mass of the halo as the companion moves away. 

6.2. Discussion 

Tidal interaction of a galaxy with its companion in a 
prograde orbit is studied using fully consistent JV-body 
models. We find that the physical properties of the 
tidally-induced nonaxisymmetric features are tightly cor¬ 


related with the interaction strength parameters P and 
S. These correlations are in qualitative agreement with 
those from the 2D results presented in Paper I, although 
the structures in the 3D models are generally weaker and 
decay more rapidly. These are also qualitatively consis¬ 
tent with the results of Elmegreen et al. (1991) who de¬ 
scribed in terms of S the formation criteria of the outer 
disk features such as double arms corresponding to spi¬ 
ral arms decoupling from the tidal tail. Byrd & Howard 
(1992) similarly found that the tidal strength P deter¬ 
mines the strength of grand-design spiral arms in inter¬ 
acting galaxies. Overall, S' is a better parameter for char¬ 
acterizing tidal tails formed in the outer regions, while 
spiral arms induced in the inner regions are better cor¬ 
related with P. 

Strength and shapes of tidal tails have been used to 
constrain the masses of dark halos (e.g., Faber & Gal¬ 
lagher 1979). Using numerical simulations, Dubinski et 
al. (1996) found that a tidal tail becomes shorter and 
less massive as the mass and size of a dark halo increase 
due to a shorter interaction time as well as a deeper po¬ 
tential well (e.g., White 1982). More quantitatively, Mo 
et al. (1998) suggested £ = (v e /v c ) 2 as an indicator of 
the susceptibility of a disk galaxy to the tail formation, 
where v e and v c refer to the escape and circular veloci¬ 
ties at R = 2 Rd, respectively. Note that £ measures the 
depth of the halo potential relative to the specific kinetic 
energy of the disk. Using galaxy models with 4 < £ < 8, 
Springel & White (1999) showed that the maximum tidal 
response T e s described in Section 3.1 is inversely propor¬ 
tional to £. They argued that tidal forcing produces a 
tail with a substantial amplitude only when £ < 6.5 in 
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Figure 17. Axially-averaged rotation velocities ( v of the primary halo (left) and the companion (right) as functions of the radial 
distance in model TA2H. The vertical arrow in each panel marks the location of the companion in the left panels and the primary halo in 
the right panels. 


their models. Dubinski et al. (1999) conducted an exten¬ 
sive survey of galaxy collisions, and confirmed that £ is 
really a governing parameter and must be less than 6.25 
for the tail formation. 

Our primary galaxy has £ = 4, corresponding to the 
lower end in the models of Springel & White (1999) (see 
their Figure 10). Based on the results of Springel & 
White (1999) and Dubinski et al. (1999), therefore, one 
can naturally expect strong tidal tails produced in all 
of our simulations. However, it turns out that a tail is 
absent in model TB3 with S = 0.07, while all the other 
models with S > 0.09 do produce a tail with strength 
depending on S. Model TB3 has a very weak tidal re¬ 
sponse at %s = 0.005 from particles captured by the 


companion, even if it is more prone, in terms of £, to 
the tail formation than any model of Springel & White 
(1999). This suggests that the tail formation should de¬ 
pend not only on £ but also on the interaction strength. 
Our numerical results imply that a criterion for produc¬ 
ing a tidal tail is S > 0.09 (or P > 0.05) for £ = 4. The 
critical values of S and P would increase with increasing 
£, since it would then become increasingly difficult to 
make disk particles climb out of the halo potential well. 

We found that spiral arms induced by a tidal encounter 
are approximately logarithmic, consistent with many of 
observed grand-design spirals including M51 (e.g., Ken- 
nicutt 1981; Elmegreen et al. 1989; Patrikeev et al. 2006; 
Shetty et al. 2007). In our models, a logarithmic shape 
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results from the fact that the spiral arms in the present 
simulations are close to kinematic density waves, slightly 
modified by self-gravity, for which the pitch angle varies 
with time and space as 


tan i = t 


-l 


d(Q — k/ 2) 
dlnP 


(16) 


(e.g., Eq. (6.26) of Binney & Tremaine 2008). The ro¬ 
tation curve of our galaxy model shown in Figure 1 
gives d(Q — n/2)/d\\\R ~ 3.5 ± 0.5 km s -1 kpc -1 at 
5 kpc < R < 15 kpc where arms are induced, making 
tan* roughly constant over R. This predicts that kine¬ 
matic density waves wind out with time as tani oc t~ l 
when self-gravity is neglected. Since f i > k/2 for galac¬ 
tic disks, this also predicts that tan i is smaller for larger 
|dS2/dlnP|, consistent with the observational results that 
the arm pitch angle has a negative correlation with the 
shear rate (e.g., Seigar et al. 2005, 2006, 2014). Our 
results show that self-gravity reduces the arm-winding 
rate to tan* a t -a75 , although it is not strong enough to 
maintain the arms quasi-steady. This holds even in 2D 
models with overestimated self-gravity at the disk mid¬ 
plane (Paper I). 

Another prediction of kinetic density waves is that the 
arm pattern speed is not strictly constant over the ra¬ 
dius but a decreasing function of R. Indeed, Westp- 
fahl (1998) found that the pattern speed of the arms 
in M81 interacting with NGC 3077 is best described 
by flp = 44.6-2.3(77/1 kpc) km s -1 kpc -1 based on the 
Tremaine & Weinberg (1984) method, after taking al¬ 
lowance for the radial variation of fl p . Meidt et al. (2013) 
applied a similar method to M51 in obvious interaction 
with NGC 5195, and found that the arms consist of mul¬ 
tiple patterns with larger pattern speeds at smaller radii 
(see also Meidt et al. 2008a,b), although it is questionable 
whether they are really multiple patterns with different 
pattern speeds or a single pattern with radially-varying 
Qp. Realistic interaction models of Salo & Laurikainen 
(2000b) for the M51/NGC 5195 pair also predict Q p close 
to the fl — k/2 curve over a range of radii where the arms 
are strong. It is interesting to note the spiral arms in 
the isolated galaxy NGC 1068 have a pattern speed that 
decreases steeply with P, indicative of a short lifetime 
of order of ~ 0.1 Gyr (Speights & Westpfahl 2011; see 
also Speights & Westpfahl 2012). All of these suggest 
that arms in real spiral galaxies are more likely tran¬ 
sient rather than possessing the characteristics of quasi¬ 
stationary density waves unless the disks are strongly 
self-gravitating. 

While we employ very simple galaxy models for tidal 
interactions, our models can well be applied to the 
M51/NGC 5195 system that have the mass ratio of 
~ 0.3-0.55 (e.g., Smith et al. 1990) and a pericenter dis¬ 
tance of ~ 20-30 kpc (e.g., Salo & Laurikainen 2000a), 
corresponding to models TAl and TA2. As Figure 7 
shows, the arm strength is maximized ~ 200 Myr after 
the pericenter passage. Lee et al. (2005) found observa- 
tionally that the age distribution of star clusters in the 
arms of M51 has a narrow peak at 4-10 Myr and a broad 
peak at 100-400 Myr. If the enhanced star formation 
is really triggered by tidal effects, the M51/NGC 5195 
pair might have undergone double interactions, the first 
one 400-500 Myr ago and the second one 50-100 Myr 


ago (Salo & Laurikainen 2000a), considering the time 
delay between the arm formation and the closest ap¬ 
proach. The peak arm strength in models TAl and TA2 
is T ~ 10 — 20%, consistent with the radially averaged 
value of ~ 15 — 20% from A-band observations (e.g., Rix 
& Rieke 1993; Rix & Zaritsky 1995; Salo & Laurikainen 
2000b; Scoville et al. 2001). The arm pitch angle at peak 
is i ~ 17° — 22° in models TAl and TA2 (Figure 10), 
which is also consistent with the observed pitch angle of 
i ~ 17.5° — 21.1° (e.g., Shetty et al. 2007; Fletcher et al. 
2011; Hu et al. 2013; Puerari et al. 2014). Based on our 
numerical results, the stellar spiral arms in M51 are in 
the nonlinear regime and their 3D density structure may 
be described by the concentrated arm model of Cox & 
Gomez (2002). 

By analyzing the properties of to = 2 spiral arms in 
the Spitzer Infrared Nearby Galaxies Survey, Kendall 
et al. (2015) very recently found that spiral morphol¬ 
ogy depends only weakly on the galaxy properties, while 
the arm strength is tightly correlated with tidal forc¬ 
ing from nearby companions. Their Figure 17 shows 
that 1Z = S TO= 2 /E m= o is an increasing function of P, 
which agrees with our numerical results qualitatively, but 
seemingly not quantitatively. That is, TZ of the arms in 
the Spitzer samples varies in the range of 0.05-0.5 for 
logP = —5-0, while our results predict 1Z ~ 0.2-1 only 
for log P ^ — 2 and no arm induced for log P ;$ — 2 (see 
also Byrd & Howard 1992). As Kendall et al. (2015) 
noted, this apparent quantitative discrepancy may arise 
because the observed P values in the Spitzer samples cor¬ 
respond to the projected distance of a nearest neighbor 
at the current epoch, while P in our models measures 
the tidal force using the 3D distance at the epoch of the 
pericenter passage. The projection effect overestimates 
the tidal forcing, while the distance at the current epoch 
is likely to underestimate the true P significantly. For 
example, a comparison between Figure 17 of Kendall et 
al. (2015) and our Figure 9 suggests that NGC 1566 with 
TZ ~ 0.2 and the current forcing estimate of log P ~ —4.3 
might probably have undergone an interaction with its 
companion NGC 1581 with logP ~ —1 at the pericenter 
passage. 

Finally, we remark on the absence of a bar in our nu¬ 
merical simulations. As Figures 5 and 6 show, tidal forc¬ 
ing in our models induces spiral structures at R xi 5 kpc, 
while the inner regions with R ^ 4 kpc remain almost 
unaffected. This is in contrast to the results of Noguchi 
(1987), Gerin et al. (1990), and Mihos et al. (1997) who 
showed that disks subject to tidal interaction can be un¬ 
stable to form a bar. The major difference between the 
galaxy models used by the present paper and their work 
is that our models possess a relatively strong bulge, while 
their models have no or very weak bulge. The strong 
bulge in our models puts an inner Lindbland resonance 
which suppresses feedback to swing amplification that 
would otherwise produce a bar near the center (Toomre 
1981). This suggests that a strong bulge may be respon¬ 
sible for the absence of a bar in non-barred grand design 
spiral galaxies. 
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